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We consider the implications of spin-phonon coupling within the slave-boson, mean- 
field treatment of the extended t-J model of a high-temperature superconductor. In 
bilayer cuprates such as YBaCuO, where the CuO-i plane is buckled, this interaction is 
linear in O displacement along the c-axis, and the coupling constant is found to be large. 
The formation of a spin singlet causes additional contributions to the phonon self-energy, 
and we calculate from these the superconductive phonon anomalies. The magnitude and 
sign of the frequency shift and linewidth broadening for various mode symmetries cor- 
respond well with Raman and infra-red light scattering experiments, and with neutron 
scattering studies. In the t-J model, spin singlet formation and superconductivity do 
not coincide in the low-doping regime, giving rise to spin-gap features and a variety of 
temperature scales in the spin response observed by NMR and neutron investigations. 
Phonon anomalies in underdoped YBaCuO compounds indeed show evidence of spin- 
gap phenomena with the same characteristic temperature, suggesting that the theory 
may offer the possibility of a unified understanding of the anomalies in magnetic and 
lattice properties. While the origin of the superconducting interaction is electronic, this 
spin-phonon coupling affords the possibility of a small isotope effect, and our estimate 
is in good agreement with recent site-selective O-substitution experiments. 
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1 Introduction 



Since the discovery of high-temperature superconductivity in oxide-based ceramics, a 
wide variety of experiments has shown strong evidence of a link between anomalies in 
the lattice structure and the onset of superconductivity in these materials. A com- 
prehensive overview of the techniques available and their results is contained in Ref. 

More specifically, a number of local probes has shown this connection: neutron 
scattering studies || |3| of the probability distribution function for atoms in the C11O2 
plane show clear local distortions at the transition temperature T c in several systems; 
ion-channeling Q shows anomalous changes in the oscillation amplitudes of Cu and O 
atoms not only at T c , but also at temperatures above it; EXAFS experiments || show 
changes in the environments of both Cu and apical O atoms; ultrasound measurements 
H show alterations of the elastic constants of the crystal. Among the more general 
probes, there have been various reports from Raman scattering and infrared reflectivity 
investigations @, || ||| which indicate anomalous temperature-dependences in the fre- 
quency shift and linewidth of phonons around and above the superconducting transition. 
A separate manifestation is the existence of an isotope effect |Tl| , albeit one which 
appears to have very different forms as a function of sample doping. Strikingly, infrared 
reflectometry measurements of both the in-plane conductivity a a ^ Jl2] and the c-axis 
conductivity a c |T3 also illustrate a connection between c-axis phonon modes and the 



properties of charge carriers. 

In this work we will be particularly interested in the issue of phonon anomalies, which 
are one direct consequence of the spin-phonon coupling under consideration. Anomalies 
in phonon mode frequencies and linewidths at T c in the Y Ba2Cu^O'j-s (YBCO) system 
have been examined since very early in the history of high-T c superconductivity [ |I4 



because of the light they can shed on both the magnitude and the symmetry (from mode 
polarization) of the electron-phonon interaction. Studies by Raman spectroscopy have 
centered largely on the 340cm -1 mode of planar oxygen which shows a particularly large 
effect in the O7 compound, and we will encounter its importance throughout the present 
work. However, it is essential to note, using the well-characterized case of this mode, 
that not all early measurements of the anomalies, performed on a variety of systems, 
reached the same result; detailed investigations of high-quality samples have since ver- 
ified the anomaly, but have shown that its magnitude is very sensitive to both oxygen 
content ||, falling sharply as the O depletion 5 increases from zero, and to impurities 
[ 15| . Because samples with non-zero S are non-stoichiometric, a considerable amount of 
disorder may be present, acting to suppress the observed effects; thus the stoichiometric 
materials Y2BCI4CU7O15 and YBa2Cu^O% will represent important benchmarks of the 
variation in anomalies with doping level. In other classes of the high-T c materials there 
remain debates about sample purity, and rather few crystals of any single compound, 
so particular care will be required in interpreting results from these. 

On the theoretical side, anomalies in the transport and magnetic properties, observed 
in both normal and superconducting states of the high-temperature superconductors 



Tq| , are considered to be manifestations of a metallic state arising near the Mott tran- 
sition due to strong correlations. While there remains no consensus on the appropriate 
theoretical description of the phenomena associated with this anomalous metallic state, 



in this work we will pursue an approach based on the t-J model [17 treated by the 
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slave-boson mean-field theory JTS, [19], |20], This framework has been shown to con- 
tain ingredients essential to an understanding of both transport properties [ 22 1 , including 
the temperature-dependence of the resistivity and Hall coefficient, and spin excitations 

, such as the different temperature-dependences 



23, 24, 25, 26, 27 
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of the 

shift and rate of nuclear magnetic resonance between high- and low-doping regions. On 



this last point, such models appear to contain a consistent description of the 
behavior |JT 



spm-gap 

noted first by Yasuoka JHJ an d subsequently the subject of much 
investigation. However, while some aspects of the link between these features and the 
phonon problem have been studied theoretically, to our knowledge there have not so far 
been any efforts to understand all of the low-lying excitations on a unified basis. 

The application of conventional electron-phonon coupling to the HTSC problem was 
made by Zeyher and Zwicknagel ||32|| , in a model which required a number of basic 
assumptions about the nature of the system, including Fermi liquid behavior of quasi- 
particles, an s-symmetric superconducting gap and a cylindrical Fermi surface. These 
authors computed the form of the superconductive phonon anomalies, and found that 
a very strong coupling constant A ~ 3 was required to reproduce the experimental 
features, which appears to be inconsistent with later studies of transport phenomena 
within the same framework 
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At the bandstructure level, the resonant frequencies 
of many phonon modes of the YBCO structure were computed in an extensive LDA 
study by Andersen et al. |34T| , who concluded that good qualitative agreement with 



experimentally observed anomalies would be given on application of the conventional 
model. These authors drew particular attention to the "dimpling" modes of planar 0(2) 
and 0(3) atoms, not only due to the strong anomaly in their out-of-phase oscillatory 
mode (340cm _1 ), but also because they cause "extended saddle points" near the Fermi 
surface of the electron dispersion in the region of the (tt, 0) points in fc-space, which 
contribute strongly to the calculations of many physical properties. Further studies also 
based on conventional electron-phonon coupling [^, |36] have given detailed descriptions 
of the Fano lineshape, and pointed out how phonon anomalies may be used to deduce 
the gap symmetry. We note briefly that in contrast to most of the conventionally-based 
treatments ||32|| , the present analysis is based on strong correlations, encoded by the 



t-J model, and takes a simple, mean-field approach containing no further assumptions 
about the nature of the system. 

We emphasize that this is not a study of a phononic mechanism for superconductivity. 
In the t-J and related models, quasiparticle pairing is the consequence of the spin 
interaction, and so its origin is entirely electronic. In this work we investigate the effects 
of phonon oscillations of the host lattice on the parameters of the model itself, showing 
that they cause an effective spin-phonon coupling. As explained in the following section, 
while the coupling constant is large, the net interaction depends also on the magnitude 
of the phonon oscillation, and is found to be small. However, its results are clearly 
observable in a variety of experimental properties, which may in turn shed much light 
on the behavior of the system itself. 

The structure of this paper is as follows. In section 2 we introduce the concept of a 
coupling between the spin and lattice degrees of freedom, and deduce the quantitative 
form of the interaction vertices. In section 3 is given a complete account of the appli- 
cation of the single-layer theory to phonon anomalies in Y BCI2CU3O7. In section 4 we 
present a qualitative discussion of the extension of the same ideas to underdoped YBCO 
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materials, leading to a connection with the mean-field phase diagram and spin-gap be- 
havior. Section 5 contains an analysis of the isotope effect, and section 6 a concluding 
discussion. Some of the key results of this paper have been presented in Ref. p7 ]. 



2 Spin-Phonon Coupling 

We consider the spin-phonon coupling arising naturally within the extended t-J model 



of a single Cu0 2 layer, which has been shown |24], g£j at the mean-field level to give a 
good account of many features of the spin excitations in both Ld2- x Sr x CuO i (LSCO) 
and YBa 2 Cu 3 7 ^s (YBCO), based on realistic Fermi surface shapes for each class of 
material. The Hamiltonian is 

H = — tijolpdjo- + JjjSj.Sj, (1) 

V> (ij) 

where the Hilbert space is that without double occupancy, corresponds to the transfer 
integrals used to reproduce the Fermi surface, and the superexchange interaction J^, 
which is usually taken to be a constant at fixed doping, is assumed to be finite only 
between nearest neighbors. 

In YBCO, the C11O2 layer is "buckled", by which is meant that the oxygen atoms 
0(2) and 0(3) lie out of the plane of the Cu atoms, as shown schematically in Fig. 1. 
Let uo + uf represent the magnitude of the oxygen displacement along the c-axis, with 
Mo the equilibrium buckling and uf the phonon coordinate, in which i refers to the Cu 
sites on the square lattice, and a to either 0(2) (a = x) or 0(3) (a = y)). Then tij 
and Jij between nearest-neighbor Cu sites have contributions linear in uf which may be 
written as 

U,i+ a = t[l-X t (uf/a)), (2) 
J i>i+a = J[l-\j(uf/a)}, (3) 

with a the distance between Cu sites, and A t and Aj the corresponding coupling con- 
stants. Microscopic estimates of A t and Aj are somewhat involved, and we proceed to 
detail the steps required in their derivation. 



The degree of buckling may be deduced from the structural data of Ref. pql , which 
shows the equilibrium O displacement to be uq = 0.256A; this data was taken at 10K, 
but the temperature-dependence is thought to be very weak, even at the superconducting 
transition. We neglect the 5% anisotropy between a- and 6-axes in the pseudo-tetragonal 
system, and thus continue also with the assumption that the atomic levels of the a- and 
7r-orbitals on 0(2) and 0(3) are the same. Although the degree of buckling is small, 
uo/a 1, its inclusion is crucial in providing a coupling which is strong and linear. 

The dependence on interatomic separation of the transfer integral t a (t n ) between 
Cu:d x 2_ y 2 and 0:p a (O'.p^) is given, following Ref. |3"5[] , by 
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Here ^Vpda (V p dn) is the transfer integral between d x 2_ y 2 (d xz ) and p a (p z ) orbitals with 



separation d = (J^ -) + u ) , directed along x according to definition, and we have 
taken the distance-dependence to be 



-7/2 



(6) 



By writing u = uq + 5u, where Su (= uf ) denotes the oscillation amplitude, for the value 
of uo above one obtains to lowest order t a (u) /t a {uo) = l—2.03Su/a and (t n (u) — t n (u )) / 
t<j(uo) = 3.065u/a. 

The final requirement is the dependence on t a and t n of Uj and Jy, and here we must 



consider in detail a model of the C11O2 plane. Following Eskes and Jefferson f40H , Jij 
can be given by the perturbative expression 



J 



4 
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2A CT + U. 



(7) 



where Ud and U p are on-site interactions for the participating Cu and O orbitals, U p d is a 
Coulomb interaction between holes on neighboring sites and Act is the charge transfer 
energy. Application of this formula alone yields a very large coupling constant A j = 10.4, 
because of the quartic power law and because the contributions from a- and 7r-hopping 
processes in the numerator combine ferromagnetically. However, use of the lowest-order 
perturbation form ([?]) is not well justified, as there are no small parameter ratios, and 
in particular the direct 0-0 hopping integral t pp will have a strong effect. A detailed 
investigation of the influence of higher-order terms 
of the C11O2 layer the effective relationship J oc tp d , 

integrals are taken to vary uniformly (t p 
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has found for the parameters 
with x ~ 2.3, when all overlap 



oc t 2/3 

V PP ^ l pd 



39fl) and A C t is held fixed. However, 



this result depends upon the variation of the transfer integrals with the bond distortion 
associated with each phonon mode, and for a mode in which t pp remains constant (A\ g 
and A 2u symmetries, section 3), the effective power is reduced to x ~ 1.7 |4T . 

The case of the hopping terms ty is complicated further by the fact that it requires 
consideration of the many inter-site transfer integrals which contribute to the hopping 
of a Zhang-Rice singlet, and not simply of the motion of a hole. This problem has been 
considered within a cell-perturbation method | j42[| , and the results for the effective power 
in t oc t pd are y 



1.0 for uniform variation, and y 



0.7 when t pp is held constant |4l 



In 



this study, our primary aim is to elucidate the effects of spin-phonon coupling in the t-J 
model at a semi-quantitative level, with attention paid to the vertex magnitudes mainly 
to ensure that the consequences discussed are not irrelevant to experiment. Thus for the 
current purposes we accept possible errors of order 20-30% in adopting the values x = 2.0 
and y = 1.0, so that Aj = 5.2 and At = 2.6 = ~Aj, and errors of similar magnitude 
incurred by neglecting any modulation of the extended singlet hopping terms. [In fact, 
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the first set are likely to lead to an overestimate, while the second, in the most naive 
approximation, would give a similar underestimate, although this depends on the mode 
symmetry.] An independent indication of the validity of these values is provided by the 
agreement they yield with the measured isotope shift in YBCO (section 5). 

In this discussion we have disregarded changes in the charge-transfer energy Act as 
the bond length varies, because we believe this to be appropriate for the deformation 
processes associated with a phonon oscillation, which are both local and screened. This 
situation is to be contrasted with that in experiments in which the bond length is made 



to vary by application of hydrostatic pressure [figl , or by atomic substitution |44], [45 
where one finds the weaker relation J oc d~ a , 4 < a < 6. In these cases long-range forces 
are brought into play, the Madelung energy of the system is altered, and it is found that 



a (^-dependence of Act is required to account for experiment . 

In the slave-boson treatment of the t-J model ([!]), the spin and charge degrees 
of freedom carried by the quasiparticles are represented by the explicit decomposition 
Q'ia — fiaty,, where fi a is a fermionic spinon, in terms of which the spin is given by 
Si = \ f\cPoLfifip-, and hi is a bosonic hole, or holon. With these operators, the coupling 
in H (P of the phonon coordinate to the spin degrees of freedom may be expressed 
in the mean-field approximation as 

-EE toOx^ ( 8 ) 

i a=x,y 

+1 JA ^(¥) [(xli+*)Xi,i+a + 2<Aj i )A iJ - + h.c] } , 

where Xi,j = E s flfjs and Ay = (f^fji - /ufa) /s/2. The first two terms, containing 
Xi,j, will contribute in the normal, or uniform RVB |[29|| , state, and will be termed 
the m-RVB vertex, while the last appears only below the singlet RVB transition, so 
constitutes an s-RVB vertex. At the onset of singlet order there will be additional 
contributions to any physical quantity due to the appearance of a finite particle-particle 
vertex represented by the 3rd term, as well as to changes in the quasiparticle propagators 
joining the particle-hole terms. In this analysis of the t-J model ]29|| we do not include a 



density-density (nifij) decoupling of the spin interaction term JSj.Sj, because it has been 
found HH to have a small coefficient in the derivation from the d-p model. In principle 



the phonon-holon coupling vertex, which is also contained in the t term, may contribute 
to physical processes, but this will be neglected because in the normal state a boson 
polarization term vanishes at q = 0, while in the Bose-condensed state corresponding to 
true superconductivity in this model the holon has no dynamics. 

While the coupling constants A t and Aj have been found to be large, the overall 
strength of the spin-phonon coupling depends also on the magnitude of the phonon oscil- 
lations u?/a, and so remains small. It can be estimated from the root mean square fluc- 



tuations of the relative magnitudes of t and J ((0) and (Q)), by using J ((U ii+a — t) 2 }/t = 
^t\J (i u ? ) 2 )' an d similarly for J. The r.m.s. fluctuation of 5u = uf is given by 
(Su 2 ) = 2Mu ' wnere is taken for the phonon frequency of interest and M is the 
mass of the O atom. Taking as a typical frequency that of the 340cm -1 mode, to 
which we will devote most attention in the following section, (5u 2 ) = (0.055A) 2 , and the 
relative fluctuations in t and J are 3.8% and 7.5% respectively. These values appear sur- 



prisingly large, but are consistent with the estimates of Haas et ah [E3], who considered 
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only the effect of a fluctuating J term on the linewidth in electronic Raman scattering. 
In what follows we will find also that the changes in physical quantities arising from the 
combined effects of these fluctuations are of the order of their square. 



3 Phonon Anomalies 

Having derived a coupling between spin and lattice modes, we consider first its modif- 
cation of phonon dynamics. In Fig. 2 are shown the four mode symmetries of a CuOi 
bilayer in which the 0(2) and 0(3) atoms may oscillate in the c-direction. We will use 
throughout the nomenclature of the symmetry obeyed by the pseudo-tetragonal 
C11O2 planar unit in the slightly anisotropic YBCO system, and refer the reader to Ref. 
14 1 for full details. Here we will be concerned mainly with the 340cm" 1 B\ g mode in 



Y Ba2Cu^O-j^ which is an out-of-phase oscillatory motion of the planar oxygen atoms 
(uf = —u![), and is also out of phase between the planes of the bilayer |14| . This mode 



has attracted experimental interest because it shows the largest effects in the supercon- 
ducting state, and there are available both detailed Raman [0, |8|, [| data, which is taken 
at wavevector q = 0, and inelastic neutron scattering results []5G|1 , which can probe all 
q. The primary quantity of experimental interest is the T-dependence of the frequency 
shift and linewidth at q — 0, in addition to which we will illustrate the dependence of the 
results on mode frequency and wavevector. The 193cm _1 B2 U mode is identical at the 
single-layer level, although in the real crystal the atomic motions are in-phase between 
layers, and so is also amenable to description by the same model at this lower frequency. 
Because this mode is infrared-silent, it may be studied only by neutron scattering, and 
it has recently been investigated in Y BCL2CU3O7 We will also show the results of 
the single-layer model for the in-phase (u* = u\ ) Ai g - and ^n-symmetric oscillations, 
and discuss the relevance of the model to experiment for these. 

3.1 Theory 

The effect of the coupling on the dynamical properties of the phonon is calculated from 
the lowest-order spinon polarization correction to the phonon self-energy H, p (q, iu n ), 
shown in Fig. 3; the full phonon propagator is given by ni^q, iuj n ) = Hq x (q, iu n ) — 
n^ p 1 (q, iou n ), where the bare propagator has the conventional form n (q, iu> n ) = 2uJo/iu 2 — 
Uq. Thus when |Ren sp | -C uj q , the superconductivity-induced phonon frequency shift is 
given to lowest order by 8uj =ReU sp , and the correction to the linewidth T is 5T =Im U sp . 
Second-order perturbation theory in terms of the interaction strength results in a fre- 
quency shift which is given at q = by 

4 ^ 1 tanh (f*T) 
6u = c(\ J J) 2 -Y / F k . ft^, (9) 

where c = (^j (5u 2 ), with (Su 2 ) defined at the end of section 2, and has the value 
c = 1.18 x 10~ 4 for the Bi g phonon. E k is defined by E k = + A|]^ 2 , with the 
spinon band energy relative to the chemical potential, and = — JA(cos k x — cos k y ) 
the singlet order parameter, which has been shown to have d x 2_ y 2 -wave symmetry in 
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the lowest-energy state in the present framework 2D| ; the negative branch given by 
— E k is illustrated in Fig. 4 for the parameters of the YBCO system used in this study. 
Finally, F k is a form factor which depends on the mode symmetry, 

F k = 2A 2 ( 7fc £ fc + x^) 2 B lg ,B 2u modes 

F k = 2A 2 V l (& + fx7fc) 2 A 2M modes, (10) 

where 7^ = cos k x + cos r/ k = cos fc x — cos k y , A = (Ay) and x = (Xy) + §7 is written 
to contain the effects of the particle-hole vertices from both J and t terms; 8 is the 
hole concentration. In Eq. ( |TUP one has i 7 ^ oc A 2 , so that there is no phonon energy 
correction due to spin coupling in the normal state: this result applies only at q = 
for the optic phonon modes under consideration. By a similar analysis, the linewidth 
broadening is given by 



8T = -c {Xjjf ^ ]T F k [8{uo - 2E k ) - 8{uo + 2E k )\ ' . (11) 

iV k 

3.2 Results 

For the following calculations we use the self-consistent solutions to the mean-field equa- 



tanh (fk) 



tions of the extended t-J model |24| for the temperature-dependence of the parameters 
(X, A, fi) for 

i) the transfer integrals appropriate to YBCO: t = 4J, t' = — |t and t" = ~t (corre- 
sponding to hopping processes between nearest, next-nearest and third- neighbors, re- 
spectively), giving the spinon dispersion shown in Fig. 4, and 

ii) a doping level 5 = 0.2, corresponding approximately to optimal doping, or Y Ba2Cu^Oi 
in this model 



We note here that in the t-J models with the superexchange J calculated in the random- 
phase approximation (RPA), the doping cannot be taken to be a parameter which may 
be varied to reproduce experimental observations. The chosen doping level is taken 
neither too near nor too far from the antiferromagnetic instability of the system, so that 
the quantities calculated may be considered to be a reasonable reflection of the effects 
of spin fluctuation enhancement, without being singular [p9[| . Thus we will not have any 
semi-quantitative statements to make with regard to doping-dependence. Throughout 
the calculations we have assumed a Lorentzian broadening of the spinon spectrum 7 = 
0.12kBT RVB , and we find that the results depend very little on value of this parameter: 
specifically, a reduction of 7 by a factor of 10 is required to enhance the frequency shift 
by a factor of 2, for mode frequencies close to the maximum shift (Fig. 5(a) below). 

^-Dependence 

Using the mean- field parameters, the frequency-dependence of 6u, given by Eq. @, has 
been evaluated numerically at T/Trvb = 0.2 and T /Trvb = 0.8, where Trvb = 0.069J 
is the onset temperature for the singlet RVB order parameter A, and the results are 
shown by the o symbols in Fig. 5 for q = 0. The change in sign with frequency for 5u 
of the Big mode is a result of the u- dependent denominator in Eq. ([|). The frequency 
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of the crossing is an approximate measure of the value of 2Afc(T) near the (tt, 0) points, 
where the gap is maximal, as these regions are favored by the B\ g symmetry. As the 
temperature is increased towards Trvb, the a;-dependence of 5ui remains qualitatively 
the same, but the sign-change occurs at lower frequencies as the gap magnitude falls. 
The experimental mode frequency u p h = 340cm _1 ~ 0.3J for the B\ g mode is near, 
but just below, the frequency of the crossing at low temperatures, so this mode can be 
expected to show a maximal effect. 

The linewidth broadening ST is shown in Fig. 6 for the same two temperatures, and 
has the form to be expected for an imaginary part of the quantities whose real part is 
shown in Fig. 5, namely a peak at the frequency of the sign-change in Su. Note that 
because the gap is (^-symmetric (point nodes on the Fermi surface) there is no region at 
low frequencies where the imaginary part vanishes. 

T-Dependence 

The temperature-dependence of 5u is shown in Fig. 7 for q = and the frequency choices 
Uq/J = 0.15 and 0.25, which are indicated by arrows in Fig. 5: the former is in the 
regime where the mode frequency is considerably less than 2Max[Afc(T = 0)] = "2A", 
and shows a sharp transition; the latter value has reasonable quantitative agreement 
with the 340cm _1 mode, and shows clearly that as the mode frequency approaches 2A, 
the shift in frequency occurs at a temperature somewhat below Trvb- This corresponds 
well to the observations of Ref. ||, where the full frequency shift develops over a 
range of temperatures below the onset. In fact, at the value u = 0.3J, closest to 
the exact experimental mode frequency u ph , there is a small, positive frequency shift 
because c^o > 2Max[A/ c (T = 0)]; however, we illustrate the anomalous effects for a mode 
frequency just below this value, as is shown to be the case in experiment, and in general 
regard the degree of correspondence between model and experimental mode and gap 
energies as a success of the initial mean-field formulation. 

Comparison with the experimental result || for the temperature-dependent broad- 
ening of the B lg mode shows extremely good agreement in functional form and in mag- 
nitude (Su ~ 5T ~ O.Olcuo) for a frequency ojq close to u p h (Fig. 8(b)). At lower 
frequencies, away from the "resonance" energy set by the low-temperature gap, ST is 
suppressed (Fig. 8(a)). Note that in Figs. 7(a) and 8(a) (and 10(a)) we retain the 
labelling Bi g for the symmetry of the phonon distortion, although in fact the frequency 
is chosen to be close to that of the B 2u mode. 

Quantitatively, the magnitude of the effects given by the model with the chosen val- 
ues of At and Aj is within a factor of 1.5 of the Raman measurements f| on YBa^Cu^O-j. 
In the light of the approximations detailed in the preceding section, such correspondence 
may be taken as justification for these, and appears eminently satisfactory within a 
mean-field treatment using no adjustable parameters. We note also that the magnitude 
of the anomalies in frequency and broadening appear to drop quickly with doping 5, and 
that neutron measurements at q = p0| suggest a somewhat smaller frequency shift 
than is seen by Raman scattering. 
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Gap Symmetry 

The frequency shift and linewidth results obtained for the B\ g and B 2u modes are con- 
sistent only with predominantly d^^a-wave singlet pairing, as this has a constructive 
combination with the phonon symmetry in the expression (|H]) for the form factor. By 
contrast, a conventional extended s-symmetric gap is found to give an immeasurably 
small frequency shift for both B- and A-type mode symmetries, because the form factors 
(Fk) are small for 5-symmetry, whil both thermal and the form factors are small for A. 
While a pairing state characterized by an intermediate phase A y /A x = e^ [31, with 



<p different from but close to n, so that the system is close to (i-symmetry, cannot be 
excluded by the present analysis alone, it appears that all forms of anisotropic s-wave 
gaps are ruled out by their antiphase combination with the phonon symmetry (see also 
Ref. H). 



q-dependence 

We have extended our analysis to the case where the phonon has a finite wavevector, and 
will consider here only the variation of the frequency shift, 5uj, as this has been measured. 
5u(q, u, T) exhibits a wealth of features as q is varied, because the quantitative effects 
are sensitive to the detailed shape of the spinon dispersion Ek (Fig. 4), as well as to the 
bare phonon frequency as illustrated above. When the phonon has finite wavevector q, 
one may show that the components of u(q) obey the relations u x (q) = ±u y (q) for every 
q as a consequence of uf = ±wf, and the expression (^) takes the more general form 

Mq) = -c (Xjf * S [f^ U (e +) - HE.)) ^ ^~ E _-l f 

- FH(i -^- /(g -» ^i%I+i)4 ' (l2) 

where the form factors for the normal and anomalous processes are given by 

F {±) = \{c\ + cf) ± ^^[( Cl A + + C 3£ + )( Cl A_ + C 3£-) 

-(c 1 e f -c 3 A + )(c 1 £--c 3 A_)], (13) 

in which the subscripts ± denote momentum labels k ± |q throughout, and (ci, c 3 ) = 
(\/2Ao7fc, xVk) for the B\ g mode, and (\Z2A r]k, Xlk) for the A 2u mode. We assume for 
simplicity that the phonon is non- dispersive, i.e. that the mode frequency changes little 
with q, and this appears to be borne out by experiment. 

In Fig. 9 are shown the frequency shifts at q = (0.2, 0)ir as a function of u at low T 
((a), cf. Fig. 5(a)), and of T at u = 0.25 J ((b), cf. Fig. 7(b)). In the cu-dependence 
there is a clear additional contribution at uo ~ 0.45J which is found to grow and then 
disappears as q x is varied between and 0.47T. This feature is readily explained by the 
fact that the strongest contributions, from the (7r, 0) and (0, n) regions, have a second 
characteristic energy separation related to the sum of the maximum gap and the energy 
of the flat part of the spinon dispersion at (tt, 0) (Fig. 4): the general denominator 
u 2 — (Ek+q/2 + Ek- q /2) 2 m ( |T2"D remains negative at these values. 
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In the T-dependence, the most distinctive feature is that the contribution to 5lu 
from the w-RVB state is no longer vanishing, so that the anomaly at Trvb-, which is 
the difference between this and the phonon self-energy in the s-RVB state, becomes 
smaller. It also appears somewhat sharper as a function of temperature, and both of 
these features agree well with the finite-g results from inelastic neutron scattering [KO] . 

In Fig. 10 is shown the difference in the self-energy correction between low (T = 
0.1T RVB ) and high (T = 1.1T RVB ) temperatures as a function of q x (a,b) and of q along 
the (1,1) direction (c). The peak features have a simple explanation on the basis of Ej~: 
viewing the integral to be performed as an exercise in maximizing the contact between 
two spinon dispersion surfaces (one inverted for the anomalous scattering contribution) , 
the initial peak comes from a considerable improvement in surface overlap around the 
(•7T, 0) points as soon as q is offset. The peak around q/ir ~ 0.3 is due to overlap 
between opposite sides of the low-energy part of the dispersion, which is spanned by 
this wavevector, or scattering processes across the "neck" of the open Fermi surface. 
Both features are strongly dependent on the exact shape of E^, and while we suspect 
that the first is unlikely to be observed in a real experiment, it is possible that the second 
will appear. The contrast between Figs. 10(a) and 10(b) suggests that the latter peak 
would only be discernible for the 340cm -1 phonon mode, where the anomalies are most 
pronounced due to the proximity of the mode frequency to "2A". Comparison with 
the experimental results of Ref . [|(| , where fewer wavevector points could be sampled, 
shows some correspondence to the decreasing trend seen in the O7 material, although 
this decrease is somewhat more rapid in the calculation. A very recent study of the 
g-dependence of 5u for the 193cm -1 B 2u mode suggests a similar g-dependence for this, 
in that it varies little with wavevector in the region of the zone center. 

3.3 Application to Other Modes 

The 193cm -1 B 2u mode (Fig. 2) is also an out-of-phase oscillation of 0(2) and 0(3) 
atoms in the plane, although in this case the motion is in-phase between the planes of the 
bilayer. This mode has recently been analyzed in YBa2Cu^Oi by Harashina et aJ. ||51|| , 
and they observe a sharp frequency shift occuring close to the superconducting transition 
temperature, with a relative magnitude 5u>/u Q ~ 1%. These features are in very good 
agreement with the results of the single-layer model shown in Fig. 7(a) for a similar 
frequency ujq. However, the authors also find that the phonon linewidth F narrows 
significantly below the transition, which leads to the model-independent conclusion of 
a strong electron-phonon coupling. Here we have taken the intrinsic linewidth to be 
a constant, independent of the spinon spectrum, because no spinon spectral weight 
is expected for optic phonon frequencies at q — 0, so the experimental result cannot 
be obtained. The correction, Fig. 8(a), from the imaginary part of the phonon self- 
energy (Fig. 3), is always positive and non-zero for a <i-symmetric gap choice in the 
approximation used. For a quasiparticle density of states corresponding to a <i-symmetric 
gap state, it is easy to argue that in the presence of a coupling, phonon modes of 
frequencies close to the maximal value of the gap will be broadened by the quasiparticles, 
while those at low frequencies will narrow, in accord with the standard picture for 
acoustic phonons ||52|| . However, in the present case a mechanism for the applicability 
of this scenario remains to be elucidated. 
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In Figs. 5-10 are shown not just results for B ig and B 2u modes, but also those for 
modes at the same frequencies with in-phase 0(2) and 0(3) oscillations, which in the 
single-layer case would correspond to Ai g (Raman active) or A 2u (infrared-active) sym- 
metry. The Ai g mode appears at 440cm -1 , above 2 A, and so if its response was similar 
to Bi g (Fig. 5(b)) would show at best a small positive frequency shift and broaden- 
ing, similar in fact to the observation of Ref. ||. The A 2u mode occurs at 307cm" 1 , 
well positioned to show strong effects. However, one sees immediately that in the cur- 
rent approximation these have negligible anomalies, a qualitative difference from the 
i?-symmetric modes which emerges from the form factors F/. (10) appropriate to d-wave 
singlet pairing. 

Experimentally, the situation surrounding the strongly infrared-active A 2u mode ap- 
pears not to be clearly understood, since neutron studies |5(| report a "mysteriously" 
small shift 8uj/ujq ~ 0.4%, whereas far infrared spectroscopy |7j] reveals a strong shift 
8uj/uq ~ 1.2%. The result from infrared reflectivity has recently been reproduced on a 
single crystal of optimally- doped YBCO ||57|| , while the neutron scattering experiment 
suffers in that this mode appears as an extremely broad peak. The failure of the present 
model to reproduce these effects is directly attributable to its single-plane nature: the 
A 2u mode develops a very strong dipole moment, so is accompanied by significant in- 
terplane charge transfer, and a bilayer formulation including interband processes in a 
weak-coupling scheme has already been shown |54| to contain this physics. Extension 
of the model to a system of two coupled planes is an area of active research, and leads 
to non-trivial questions about the allowable symmetries of the superconducting gap in 
a bilayer which are themselves of intrinsic interest. In the case of the A\ g mode, its 
atomic motions will also cause charge motion within the Cu0 2 plane, the treatment of 
which is beyond the basic t-J framework. 

Finally, the effects studied will be strongly suppressed in the ^-symmetric phonon 
modes of 0(2) and 0(3), where atomic displacements are parallel to the plane, as the 
modulation of Uj and Jy will be only quadratic in the phonon coordinate (c/. (@,|3D). In 
addition, we would not expect to find significant superconductive anomalies in phonon 
modes involving motions of atoms in the unit cell which are not located in the Cu0 2 
layers, with the possible exception of apical 0(4), which is strongly coupled to the planar 
system as indicated in Ref. E|. In large part these qualitative expectations are borne 



out by experiment, where no other strong anomalies are observed, except indeed for a 
moderate feature in the 500cm _1 ^-symmetric mode of 0(4) ||, and a very curious, 
low-frequency mode in some YBCO compounds which appears to involve the c-axis 
motion of the Ba atom (which we note lies also in the plane of apical O). 

A concluding comment is in order on the possibilities for phonon anomalies in other 
classes of high-T c ceramics. The most extensively-studied groups of compounds are 
the La 2 - x Sr x CuO/i and related Nd 2 _ x Ce x CuO<± series with variable doping x, and the 
parent R 2 CuO^ series with a selection of elements R; in these cases there is now well- 
documented evidence for the link between superconductivity and local structural anoma- 



lies [|5£| . However, neither spectroscopic nor neutron investigations have revealed, in the 
modes examined so far, any clear signature of phonon anomalies of the type we 
consider, even in the 5-symmetric distortion of the single Cu0 2 layer in the unit cell 
(although this occurs above the characteristic value of "2A"). While this result is exactly 
that predicted by the theory, because the Cu0 2 plane is essentially unbuckled above T c 
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(HTT phase), so that phonon modulation of the t and J terms would be quadratic, and 
thus small, it is by no means certain that the observations correspond to an intrinsic 
property of the materials, rather than being an effect of impurities (section 1) or, more 
fundamentally for Raman studies, the lack of true inversion symmetry. In the LTO and 
LTT phases, the plane does undergo a distortion, consisting of alternating tilts of the 
CuOq octahedral units in the [1,1] and [1,0] directions, respectively, but against this 
background equilibrium configuration of positive and negative u values (section 2), the 
net coupling for any combination of A- and 5-symmetric gap and phonon symmetries 
at q = will cancel at linear order. 

The BSCCO group of compounds also shows a buckling distortion. While exact 
structural determinations are complicated by the presence of dislocations, modulations 



and disorder, there have been satisfactory characterizations of Bi^Sr^CaCu^O^ [BO, 61 



which show the nature of the CuO<i plane buckling to be that of the LTO distortion of 
the single plane in LSCO compounds. From the theory, we would expect no observable 
anomalies due to the cancellation of linear terms in such a structure (above), and this is 
consistent with experimental reports p2|, which indicate that a frequency shift for the 



336cm -1 Big mode is absent in Bi2Sr2CaCu20s-rela,ted materials. Given the sensitivity 
to impurities of the anomalous effects documented in YBCO materials, these results are 
yet to be verified. There remains no data on phonon anomalies also for the Tl- and 
ifg-containing materials, although there has been a detailed investigation of the local 
structural anomaly in the former H. 



4 Underdoped Phase and Spin Gap 

In the preceding sections we have detailed the predictions of a model where the phonon 
anomalies are coupled to spin singlet formation, which gives as their onset temperature 
not the superconducting critical temperature T c but the s-RVB condensation tempera- 
ture Trvb- This result raises the interesting possibility of probing the spin-gap behavior 
found in members of the YBCO class in the low-doping regime, by which is meant dop- 
ing levels below that required for the optimal T c , by considering the frequency shifts of 
particular phonons. The definition of the spin gap varies among authors and experi- 
ments, but we take it to mean the loss of spectral weight in the spin response which sets 
in at some temperature Tq above T c , and has been best characterized by observations of 



the temperature-dependence of the NMR relaxation rate [p0[ , which is maximal at this 

There exist already several experimental reports |7], |8|, |9| of anomalies in the frequency 
shift well above T c , whose onset temperature corresponds closely to that where the 
NMR rate exhibits a maximum. Recent, highly accurate Raman scattering studies of 
phonon anomalies in the stoichiometric underdoped YBCO compounds 



and Y BaiCu/±0% |SH| add considerable weight to these considerations as they show 
clearly the onset of a frequency shift at some T ~ 150K, followed by growth of this 
shift as temperature is lowered, until a saturation below T c , as shown schematically in 
Fig. 11. Such features may be understood on the basis of the mean-field phase diagram of 
the extended t-J model f29|, reproduced in Fig. 12, in which Trvb is indeed higher than 



T c only in the low-doping region. In this framework, there is an onset of short-ranged 
singlet RVB order around T , which is then identified with the crossover temperature 
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Tjivb [ S2|| ) and an increasing correlation length of the coherent regions with decreasing 



temperature until fully coherent, long-range order is attained at T c . The saturation of 
the magnitude of the phonon frequency shift below T c is a consequence of the finite- 
frequency nature of the excitation (Fig. 7(b)). We speculate that studies of phonon 
anomalies at higher wavevectors might show full development of the energy correction 
at temperatures closer to T , but on the basis of Fig. 10 doubt that these would be 
measurable. 

These results for the change in nature of the transition between low- and high- 
doping regimes are consistent with the observation that the specific heat anomalies at 
T c in each case are qualitatively different p3l B4|. The existence of two temperature 



scales may provide an explanation for the contrasting low-g behavior of the energy 



shifts in nominal 06.92 and O7 compounds observed in Ref. |p0[ : here the authors 
took as their high temperature for comparison a value of 1007^, which if the lower- 
doped compound were to possess a spin gap above T c could very well remain in the 
regime with a phonon anomaly, i.e. below T . We note also that there have been some 
experimental reports which indicate a true transition, as opposed to a crossover, around 
the higher temperature in several material classes, notably low-doped phases of YBCO, 
and it is possible that this may be identified with a lattice instability |6|, |65|, |66fl . 

In conclusion to this section we comment that the present theory contains only one 
temperature scale, Trvb, and thus a detailed quantitative explanation of the low-doping 
regime remains beyond its scope (see the related note in the introduction to section 
3). Within the same framework, a more accurate account of the features reproduced 
in outline here will require additional physics, perhaps in the form of an improved 
treatment of the boson degrees of freedom, to restore the lower temperature scale of the 
true superconducting transition. 



5 Isotope Effect 

The isotope effect was one of the critical pieces of evidence which pointed the way to 
the formulation of the BCS theory for superconductivity in conventional metals. In the 
high T c materials the situation is somewhat complex, and a comprehensive review with a 
wealth of experimental evidence and a summary of theoretical models to date is given in 
Ref. |J7| . In brief, the qualitative trend in each material class appears to be a situation 



where the oxygen isotope effect, parameterized by a = —d\nT c /dln.M ~ — ^^^7, 
is close to vanishing at the optimal doping level, but rises smoothly to values of at 
least a = 0.5 on sufficient under- and over-doping. The question of sample purity and 
randomness away from the stoichiometric compounds (section 1) may be pertinent here 
once again, but this form has been observed in almost all types of system. In general, 
electronic models such as the current one do not address the question of the isotope effect 
at all, and in this work we will have the modest aim only of elucidating the influence of 
the spin-phonon coupling at the near-optimal doping level. 

The thermodynamic properties of the coupled system may be studied self-consistently 
by including phonon-spinon terms in the free energy F (Fig. 13(a)). In a Ginzburg- 
Landau expansion, the free energy difference between a normal state and a superconduct- 
ing one with a single order parameter A has the form Fs — F^ = A(T)A 2 + B(T)A A + . . . 
It is tempting to postulate that the A 2 contribution of the additional parts acts to en- 
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hance the stability of the spin singlet state, providing a natural linkage of the spin gap 
to the lattice structure which might account for the variety of characteristic tempera- 
tures discussed in the previous section. However, the relative magnitude of the effective 
phonon-spinon vertex may be estimated as shown at the end of section 2, and because 
the contributions due to fluctuations in t and J require two such vertices, their effect 
can be expected to be of order 10~ 2 . Here we estimate this effect by evaluating the A 2 
contributions from the lowest-order phonon-spinon coupling term at T^vb, using the 
spinon propagators for the normal, or w-RVB, state. 

Omitting details of the lengthy but straightforward calculation, the free energy con- 
tributions from phonon coupling are shown in Fig. 13(b) for the terms (i), the s-RVB 
vertex part, (ii), a pair of diagrams corresponding to a spinon self-energy correction by 
the phonon, and (iii) phonon vertex correction of the u-RVB vertex. Near the transition, 
Fs~Fn may be linearized in temperature, so that the coefficient of the A 2 term becomes 
A(T) = a(T — Trvb)- The phonon contributions are approximately T- independent, so 
have the schematic form — 6A 2 , and can be written as a constant shift ST^vb — b/a; 
in Table 1 are given the values of the relative shifts STrvb/Trvb arising from each of 
the diagram types in Fig. 13(b), and for the four phonon symmetries and frequencies 
(Fig. 2). One observes from the figures that the anomalous vertex part remains ap- 
proximately constant for all of the modes, while each of the other two parts is small 
for the A-symmetric modes, but has large contributions from the 5-symmetric ones, as 
would be expected on the basis of the phonon anomalies computed in section 3. The 
figures can be shown to be approximately constant by performing the same calculation 
at different temperatures, using the slightly different self-consistent values of the u-RVB 
order parameter \ an d chemical potential \i at each. The total free energy shift for these 
four modes of a bilayer is compared with twice the free energy of the single-layer system 



29| . The net contribution corresponds to a 1.1% enhancement of T RV b, a magnitude 
exactly in line with the qualitative expectation, which to the same degree of accuracy is 
a 1.2K enhancement. The isotope shift for 16 — > 18 O on this quantity is approximately 
-6%, leading to the prediction of a total shift due to in-plane oxygen of -0.07K. The 
equivalent a parameter is 0.005. We would expect the quantitative agreement to be 
improved when interplane hopping processes are taken into account (section 3), in order 
to describe correctly the contributions from the A-symmetric modes (Table 1). 

This result may be compared with a recent site-selective substitution study of YBCO 
near optimal doping, performed by Zech et al. [171]. These authors were able to achieve 



high exchange rates of 16 by 18 during crystal preparation, while discerning with good 
resolution the location of the exchange, whether it was on planar 0(2) and 0(3) sites, 
or apical 0(4) and chain 0(1) sites. They measured a total O isotope shift (complete 
exchange) of -0.25K, of which at least 80%, or -0.20K, was found to be due to modes of 
in-plane oxygen. Qualitatively, this conclusion contains important agreement with the 
current picture, that it is the planar O modes which are most strongly coupled to the 
superconducting interaction. At the quantitative level, the theoretical result from the c- 
axis modes appears to be too small by a factor greater than 2, which, while constituting 
an acceptable level of agreement given the nature of the model, is not in as good accord 
as the phonon anomaly results within the same approximation. 

We have shown that a realistic spin-phonon coupling within an electronic model 
for superconductivity can show the small isotope effect consistent with experimental 
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observation in the high-T c superconducting ceramics at optimal doping. In the present 
theory we would not expect to see significant alterations in the magnitude of this shift 
on changes in doping, which will act to move the chemical potential and cause minor 
changes in the spinon dispersion (Fig. 4) away from the optimally-doped level. From the 
phase diagram discussed in section 4, one observes that T c falls with increasing doping 
in the overdoped region, where it is given by Trvb, and with decreasing doping in the 
underdoped region, where it is given by T B (Fig. 12). With AT C (AT RVB ) approximately 
constant, the threefold reductions in T c achieved experimentally by overdoping in some 
systems [67] would give a commensurate increase in a; in the underdoped regime, the 
changes in T# with doping have not been estimated. 

From these considerations, it seems unlikely that the t-J framework alone contains 
the origin of the spectacular changes in isotope effect with doping reported in experi- 
ment, or that a purely electronic picture could account for the apparent changeover to 
large values (~ 0.5) of a which these data present. However, the situation is complicated 
due to questions of O disorder, and even equilibrium position, on doping and isotope 
substitution, the actual fraction of O in the plane sites, and the sensitivity of charge 
fluctuations to these uncertainties. At minimum, we may suggest that the results show 
once again the spin-phonon interaction to be a detectable effect which can be used to 
probe the physics of the high-T c materials, but to be only a small perturbation on the 
dominant processes. In the case of the isotope effect, this perturbation is seen clearly at 
the optimal doping level where the dominant contribution vanishes, but a full explana- 
tion of the latter will require additional or separate physics. Ideas which have been put 
forward |67] include van Hove singularities, interlayer pair tunneling, bipolaron forma- 
tion, strong Coulomb correlations and isotope-dependent hole concentration, but many 
models based on these contain a conventional, often phenomenological, electron-phonon 
interaction, and suffer from other weaknesses which this entails. 



6 Conclusion and Discussion 

We have proposed a theory of spin-phonon coupling, based on the mean-field approxima- 
tion to the extended t-J model of a single CuOi plane, which provides good agreement 
with a variety of experimental observations on high-T c cuprates of the YBCO class. A 
key element is the buckling of the C11O2 plane present in these materials, which causes 
the coupling to be linear in O displacement along the c-axis. The model contains the 
following features: 

i) it gives, without recourse to parameter-fitting, a semi- quantitative account of the 
frequency and linewidth anomalies observed in the important and well-characterized 
Bi g - and I^-symmetric phonon modes of planar oxygen atoms in the stoichiometric O7 
compound. 

ii) the symmetry of the superconducting gap consistent with experiment is found to be 
predominantly d x 2_ y 2. 

iii) the conclusions are supported by the results obtained for different phonon mode 
types, and from observations in other materials. 

iv) those features of the phonon anomalies which are not given by the present treatment 
are expected to be readily described within the same framework by consideration of a 
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bilayer structural unit, and by a more sophisticated treatment of the intrinsic phonon 
linewidth. 

v) at lower doping levels, away from that for optimal T c , there is good qualitative 
agreement, which provides additional strong evidence, consistent with the results from 
NMR and inelastic neutron scattering, for the origin of the spin gap, 

vi) there is a satisfactory degree of correspondence with the measured isotope shift near 
optimal doping, and the phonon modes giving the most important contributions to this 
are identified. 

The single-layer t-J model for the strongly-correlated CuOi system has thus been shown 
to contain an additional class of physical phenomena, namely those related to structural 
anomalies. The theory is found to be well-suited to illustrating a variety of interesting 
properties observed in experiment, and to constitute a basis for further development. 

Based on the results for alteration of phonon dynamics by spin excitations, it would 
seem logical also to seek the alteration of spin dynamics by their interaction with the 
phonon degrees of freedom. In fact this study was motivated initially by the failure 
of the extended t-J model to reproduce at the mean- field level some features of the 



YBCO spin excitation spectrum, notably the AlmeV peak in the O7 compound ftq, 59 



The close coincidence of this resonance energy with the frequency of the B\ g mode, 
whose superconducting anomalies are the strongest, and the T-dependence of the peak 
intensity, which scales with A 2 oc p s , the superfluid density, make the present spin- 
phonon coupling theory appear well-suited to explain these features. However, as in 
sections 3 and 5, the phonon contributions to the spin susceptibility through the spin- 
phonon coupling vertices of section 2 are of order a few percent, while the AlmeV peak 



dominates the entire spin response [70]. One may consider instead direct coupling of 



the phonon to spin-flip processes via the spin-orbit (Dzyaloshinskii-Moriya) interaction 



in the buckled CuO<i plane |7I[, but this will be 2 orders of magnitude smaller still. We 
note also that the exact position of the resonance depends on the doping 5, with the 
energy declining on moving to lower doping |7"2 |, while the frequency of the B\ g mode 
appears to remain approximately constant. In addition, there is no trace of an analogous 
feature in the spin response associated with the B 2u phonon mode, which also shows a 
considerable anomaly. 

While we have explored in some detail the consequences of a spin-phonon coupling 
in this class of model, many questions remain. One of the major shortcomings of the 
theory is its failure to deal systematically with a physical range of hole-doping levels, as 
explained in sections 3 and 5, but despite this it does provide a consistent description of 
the spin gap, as discussed in section 4. Another issue which remains a drawback of the 
slave-boson formulation is the treatment of the holon degrees of freedom: in the model 
they have an intrinsic tendency towards Bose condensation, and thus contribute little to 
the dynamical phenomena, even within a gauge-field scheme where the appropriate spin- 
and hole-backflow requirements are enforced, leading to the Ioffe-Larkin compostion 
rules for transport properties [731 Ell . Some studies indicate that the holes also possess 



a Fermi surface |74], |75[, which suggests that the appropriate physical description will 
require hard-core bosons or secondary statistical transmutation. In the present context, 
this would lead to the introduction of a meaningful second temperature scale in the 
low-doping regime (section 4), and the possibility of a realistic account not only of the 
development of the phonon anomalies, but also of the features in the spin response 
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detailed by NMR and inelastic neutron scattering experiments. 

To elaborate on the observable consequences of a gauge-theory approach to describing 
fluctuations around a mean-field solution, we have considered the possibility of gauge- 
field screening [22] of the phonon anomalies. This is discussed briefly in the appendix 
for the massive, longitudinal gauge modes, which can in principle have a Fano-type 
coupling to phonons. We show that the corrections due to gauge fluctuations are small or 
vanish, depending on the phonon mode symmetry. However, the effects of the massless, 
transverse gauge modes in providing quasiparticle self-energy and vertex corrections have 
been studied only briefly in the context of their influence on the dynamical properties 
of the spin degrees of freedom ||76|| , and remain an open question. 



In conclusion, we propose that the model discussed represents a useful step in the 
formulation of a coherent theory of spin excitations, transport properties and lattice 
degrees of freedom. The agreement between theory and experiment achieved at the 
current level is encouraging in efforts to construct a unified picture of the anomalous 
metallic state, whose characterization is vital to a full understanding of high-temperature 
sup er conduct i vity. 
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Appendix 

Fluctuations about the mean-field solution of the order parameters in a model of the type 



we consider may be taken into account by the introduction of a gauge field |22|, whose 
transverse components, corresponding to phase fluctuations, turn out to be massless, 
and whose longitudinal component, corresponding to density fluctuations, is massive. 
Without going into the full details of the gauge-field formulation, we may consider the 
effect of possible gauge-field screening in the spin-phonon coupling problem as follows. 

The bare gauge field is merely a restatement of boson-fermion coupling within the 
slave boson decomposition, and has no dynamics, but the boson and fermion polariza- 
tion terms, occurring in sequences of any length, act to generate effective dynamical 
properties p2| . Restricting the discussion to Fano-type processes, the full gauge- field 
correction by polarization terms to the lowest-order spinon-phonon diagram shown in 
Fig. 3 is then simply the diagram in Fig. Al(a), where the dotted line is the full dynam- 
ical gauge field propagator — l/(n B + lip). Because the phonon vertex, from t and J 
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terms, contains only spin-spin and density-density components, it may couple only to the 
longitudinal part of the gauge field, whose propagator is simply that of a mass M ~ J, 
while the magnitude of the gauge-field vertex is approximately J; these considerations 
may be made more quantitative, but the important feature of the result is unaffected. 
As only the spinon coupling to the phonon is considered, each end "bubble" in the full 
gauge-field-mediated sequence, with one fc-dependent phonon and one constant gauge 
vertex, may be denoted by II' , in which case the total corrected expression at the order 



Brief examination of the diagram (Fig. Al(b)) for U' sp shows that at q = the 
expression must be 0(A 2 ): with a w-RVB vertex the anomalous propagators are re- 
quired, while with an s-RVB vertex, of 0(A), a single anomalous propagator appears, 
so the correction will unavoidably be 0(A 4 ), and unimportant. One may write down 
an expression analogous to @ for Re IT ' and find that the form factor is 



This is odd in rj^ for the .B-symmetric modes (in fact for both s- and c/-symmetry of 
the singlet order parameter), and so vanishes upon /c-integration, while for A-symmetric 
modes it remains of the order of n sp , and so becomes irrelevant on taking the square. 
Thus the consequence of considering fluctuations through a gauge-field coupling, which 
will involve only massive longitudinal modes with no singularities in momentum space, 
is that their correction is small or vanishing. 



of Fig. 3 is Tl sp - IT 




B lgj B 2u modes 
Ai g , A 2u modes. 



(Al) 
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Figure Captions 



Fig. 1: Schematic representation of the Cu0 2 layer, a) Planar structure, b) "Buckling" 
deformation of the equilibrium positions of 0(2) and 0(3) atoms out of the plane of 
the Cu atoms, appropriate for YBCO. a, b and c represent the crystal axes (assuming 
tetragonal symmetry). 

Fig. 2: Schematic representations of phonon modes involving c-axis oscillations of in- 
plane oxygen atoms in the YBCO bilayer structural unit: a) Ai g , b) Bi g , c) A 2u and d) 

B 2u . 

Fig. 3: Diagrammatic form of the lowest-order contribution to the phonon self-energy 
U sp due to coupling to spinons. The thick line denotes the spinon propagator in Nambu 
representation, and the dot the corresponding 2x2 spinon-phonon vertex. 



Fig. 4: Energy dispersion — — + for a YBCO-like system, where the transfer 
integrals in the extended t-J model are t = 4 J, t' — — ~t and t" = |t, corresponding 
respectively to nearest-, second- and third- neighbor hopping processes. 

Fig. 5: Phonon frequency shift 5uj for Bi g (o) and A 2u (x) modes at q = as a 
function of frequency at a) T = Q.2Trvb and b) T = O.STrvb- The arrows indicate the 
frequencies whose T-dependence is illustrated in Fig. 7. 

Fig. 6: Phonon linewidth broadening ST for B\ g (o) and A 2u (x) modes at q = as a 
function of frequency at a) T = 0.2Trvb and b) T = 0.8Trvb- The arrows indicate the 
frequencies whose T-dependence is illustrated in Fig. 8. 

Fig. 7: Phonon frequency shift 5uj for B\ g (o) and A 2u (x) modes at q = as a function 
of T for mode frequencies a) u = 0.15 J and b) u = 0.25 J. 

Fig. 8: Phonon linewidth broadening 6T for B\ g (o) and A 2u (x) modes at q = as a 
function of T, for mode frequencies a) u> = 0.15 J and b) u = 0.25 J. 

Fig. 9: Phonon frequency shifts Sou for B\ g (o) and A 2u (x) modes at wavevector at 
q = (0.2,0)7r; a) as a function of frequency at T = 0.2Trvb and b) as a function of T 
for mode frequency ojq = 0.25J. 

Fig. 10: Phonon frequency shift Sou for Bi g (o) and A 2u (x) modes of wavevectors a) 
q x in the (1, 0) direction, for a phonon of bare frequency u = 0.15 J, b) q x in the (1, 0) 
direction, for a phonon of bare frequency cuq = 0.25 J, and c) qj\[2 in the (1, 1) direction, 
for a phonon of bare frequency u = 0.25 J. 

Fig. 11: Schematic representation of the frequency shift of the 340cm -1 B\ g phonon 
mode as a function of temperature in the underdoped compound YB^Cu^Os, from the 



data of Ref. |pb| . T c = 82i^ is the temperature of the superconducting transition, and 
T sg ~ 1507^ denotes the onset of the anomalous frequency shift. 

Fig. 12: Mean-field phase diagram for the extended t-J model in the slave-boson formu- 
lation, showing the existence of a spin gap at doping levels below that for the optimal 
T c (after Ref. 
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Fig. 13: Diagrammatic representation of a) the free energy contribution due to the 
lowest-order spinon-phonon couping term, and b) those parts of order A 2 , to be evalu- 
ated in the normal state at Trvb- 

Fig. Al: Gauge-field correction to the lowest-order spinon polarization, a) Represen- 
tation of all orders of correction due to gauge field, b) End bubbles in diagrammatic 
sequence, showing 0(A 2 ) nature. 
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(i) 


(ii) 


(iii) 


total 




1.545 x 1(T 3 


3.453 x 1(T 4 


-3.588 x 1(T 4 


1.532 x 1(T 3 


Big 


1.319 x 10" 3 


-1.775 x 10" 4 


3.210 x 10" 3 


4.351 x 10" 3 


A 2u 


1.812 x 10" 3 


3.755 x 10" 4 


-4.223 x 10" 4 


1.765 x 10" 3 


B2u 


1.718 x 10" 3 


-2.558 x 10~ 3 


4.547 x 10" 3 


3.708 x 10" 3 



Table 1: Lowest-order spinon-phonon coupling contributions to the s-RVB transition 
temperature. The magnitude is relative to Trvb = 0.069 J, and the sign is chosen so 
that a positive contribution corresponds to an enhancement of Trvb- Figures are given 
at T RV b for the 4 modes involving c-axis motion of planar 0(2) and 0(3) atoms (Fig. 2), 
and for the 3 types of process shown in Fig. 13(b). 
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